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1.  INTRODUCTION  AND  SUMMARY 


By  dimensional  analysis,  one  can  establish  that  gravity  plays  an  increas- 
ingly important  role  in  cratering  and  ground  motion  for  larger  yield  events,  and 
that  the  effects  of  material  strength  are  relatively  less  important  at  larger  yields. 

It  has  long  been  recognized^  that  these  considerations  prevent  sealing  high 

explosive  phenomena  up  to  multi-kiloton  yields  and  higher  bv  the  so-called 

1/3 

hydrodynamic  scaling  (W  ). 

The  data  base  confirms  the  inapplicability  of  scaling  crater  dimensions 

1/3  (2-111 

by  W .A  considerable  amount  of  research  has  been  devoted  to  develop 

scaling  laws  using  other  yield  exponents.  Many  of  these  research  efforts  have 

specifically  recognized  the  influence  of  gravity.  Scaling  by  means  of  a yield 

(12) 

exponent  remains  controversial  , however,  and  it  is  highly  desirable  to  be  able 
to  simulate  large  yield  events  in  various  media  by  an  inexpensive  means.  As 
discussed  in  Section  2,  the  effects  of  gravity  may  be  taken  into  consideration  by 
conducting  the  experiment  in  a centrifuge.  If  the  linear  dimensions  of  an 
experiment  are  scaled  down  by  the  same  factor  by  which  the  gravitational 
acceleration  is  scaled  up,  dynamic  similitude  is  preserved.  That  is,  if  one  uses 
the  same  material  in  both  experiments,  its  material  properties  scale  correctly  as  < 

long  as  rate  effects  are  unimportant.  This  comment  applies  to  the  material  « 

properties  of  the  geologic  medium,  the  high  explosive  and  air.  Thus,  both  the 
direct-induced  and  the  airblast-induced  ground  motions  scale  properly  in  a centri- 
fuge experiment.  It  can  be  shown  that  the  yield  scales  as  the  cube  of  the  1 

experimental  dimension,  so  that  at  500  g's  one  gram  of  HE  simulates  the  effects 
of  125  tons  of  the  same  HE. 

i The  work  described  in  this  report  was  undertaken  to  provide  theoretical 

support  and  analysis  for  an  experimental  program  conducted  by  R.  M.  Schmidt 
of  the  Boeing  Company  under  separate  contract  to  the  Defense  Nuclear  Agency. 

We  should  point  out  that  the  final  analysis  of  these  test  results  will  be  given  by 
Schmidt  in  a separate  report,  and  the  results  as  discussed  here  are  preliminary. 

It  is  not  expected  that  our  conclusions  will  be  modified  by  additional  data 
reduction. 
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A brief  review  of  previous  applications  of  centrifuge  testing  to  cratering 
and  related  problems  is  given  in  Section  2.  The  development  of  a high  explosive 
charge  equivalence  for  the  JOHN1E  BOY  event  is  discussed  in  Section  3.  This 
charge  was  incorporated  into  the  Boeing  test  matrix.  Finally  in  Section  4,  the 
results  of  the  Boeing  tests  are  discussed  within  the  context  of  the  results  of 
other  cratering  events.  The  program  objectives  and  the  conclusions  and  recom- 
mendations which  we  feel  have  resulted  from  this  effort  are  presented  in  the 
remainder  of  this  section. 

1.1  PROGRAM  OBJECTIVES 

The  ultimate  goal  of  centrifuge  simulation  of  large  yield  events  is  the 

reproduction  of  the  crater  results  observed  on  the  nuclear  events  at  the  Pacific 
(13) 

Proving  Ground  , and  hopefully  some  experimental  guidance  for  resolution  of 
the  differences^1^  between  empirical  predictions  and  detailed  finite  difference 
calculations  of  craters  from  nuclear  surface  bursts.  Funding  limitations  imposed 
more  modest  objectives  on  this  first  effort.  The  Boeing  test  matrix  was  limited 
to  twelve  events.  The  first  series  of  six  shots  was  conducted  in  Ottawa  sand 
so  that  extensive  testing  performed  by  Piekutowski^1^  could  be  used  as  a tie 
point  to  results  at  normal  gravity.  The  objectives  of  these  first  six  tests  were: 

• Determine  reproducibility 

* 

• Quantify  centrifuge  scaling 

• Explore  effects  of  particle  size  variation. 

A Program  Review  Meeting  was  held  at  the  end  of  this  phase  of  testing. 
It  was  the  consensus  of  the  group**  that  while  the  objectives  of  this  first  series 
* of  shots  had  been  achieved  quite  well,  Ottawa  sand  did  not  represent  a satisfactory 

material  for  simulation  of  events  in  desert  alluvium.  Field  material,  with  the 
larger  particles  removed,  was  provided  to  Boeing  by  R.  W.  Henny  of  the  Air 
Force  Weapons  Laboratory.  The  last  six  tests  were  conducted  in  this  medium 
and  devoted  to  the  following  objectives: 

• Verify  reproducibility  and  quantify  centrifuge  scaling 
in  this  new  material 


♦See  Section  2.2.3 

♦♦Included  representatives  of  DN A/SPSS,  AFWL,  PacTech  and  Boeing. 
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• Test  the  predicted  JOHNIE  BOY  simulation 
charge 

• Explore  the  sensitivity  of  the  simulation 
charge  to  depth  of  burial. 

1.2  CONCLUSIONS 

The  finite  difference  calculations  of  the  early  ground  motion  for  the 
JOHNIE  BOY  event  indicated  that  the  best;  small-scale  simulation  with  PETN 
would  be  a sphere  with  a depth  of  burial  of  0.348  cm  of  a 1.2  gm  charge 
conducted  at  345  g.  The  results  of  this  simulation  were  low  on  crater  volume 
by  22%  and  high  on  aspect  ratio  by  20%.  This  variation  is  consistent  with  the 
uncertainty  in  the  details  of  the  JOHNIE  BOY  emplacement'  , the  material 
properties  at  the  time  of  the  event,  and  the  possible  variation  in  the  detonation 
characteristics  of  the  PETN  charge.  Increasing  depth  of  burial  by  approximately 
one  charge  radius  led  to  a crater  volume  that  was  high  by  16%  and  slightly 
reduced  the  aspect  ratio. 

The  data  from  the  Boeing  tests  in  both  materials  (Ottawa  sand  and 
desert  alluvium)  confirm  the  exchangeability  of  the  charge  mass  and  the  cube  of 
gravity,  one  aspect  of  centrifuge  scaling  (i.e.,  a lack  of  rate-dependent  effects), 
and  no  significant  perturbations  in  crater  shape  due  to  Coriolis  accelerations  were 
observed.  However,  it  was  found  that  crater  volume  in  Ottawa  sand  was 
proportional  to  rather  than  to  Y.  The  tests  in  alluvium  indicate  that 

apparently  small  changes  in  test  bed  density  and  moisture  content  can  lead  to 
substantial  variation  in  crater  dimensions.  While  these  changes  are  large  enough 
to  be  observed  in  well-controlled  laboratory  experiments,  they  represent  the 
scatter  likely  to  occur  in  field  events.  By  understanding  these  variations,  the 
proper  bounds  can  be  established  above  and  below  "nominal"  cratering  curves  used 
for  targeting  and  vulnerability  analyses. 

Based  on  analysis  of  crater  dimensions  from  the  centrifuge  tests,  and 
on  comparison  with  other  near-surface  bursts,  we  conclude  that  several  empirical 

relations  can  be  written  expressing  systematic  variations  in  crater  dimensions. 

4 5 3 

For  scaled  yields  (see  Eq.  (19),  Section  4.2)  less  than  about  10  or  10  g’s  kg, 
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zero  height  of  burst  craters  in  Ottawa  sand  and  in  reconstituted  alluvium  fall  on 
separate  curves.  For  non-cohesive  dry  material  (e.g.,  Ottawa  sand)  crater  volume 
V will  be  within  ten  percent  of 


(2.66  x )0'4m3)  -0.477  w0.841 

P 9 


(U6S), 

(U66), 


-0.159 


3 2 

where  P is  the  density  in  kg/m  , g is  the  acceleration  in  g’s  (9.8  m/sec  ),  W is 

the  charge  mass  in  kg,  and  U and  6 are  the  Chapman-Jouget  (C-J)  particle 

velocity  and  density,  respectively,  with  the  subscripts  indicating  C-J  properties 

for  TNT  and  the  explosive  used.  For  moderately  cohesive  soils  (e.g.,  reconstituted 

KAFB  alluvium  and  natural  soils  for  many  high  explosive  field  tests)  crater 

volumes  obey 


( 6 . 2 +-  2.2)  x 10~5m3  -0.070  w0.977 

P 9 


(U66) 
_ J 

(U6<5 ) 


-0.023 


There  are  no  data  for  high  explosive  craters  in  alluvium  at  scaled  yields  higher 
5 3 

than  about  10  g's  kg,  but  it  is  expected  that  they  would  have  volumes  closer 

to  (1)  than  to  (2).  Craters  at  non-zero  heights  of  burst  will  have  volumes  given 

0.279 

by  (1)  or  (2)  multiplied  by  10  0-115H/W  where  H is  the  height  of  burst  in 

0 28 

m.  If  the  same  W ' height  of  burst  scaling  is  applied  to  nuclear  craters,  it 
is  found  that  crater  volumes  from  lkt  nuclear  bursts  in  alluvium  and  all  of  the 
Pacific  nuclear  craters  fit  a single  height  of  burst  and  yield  relation 


V = V y°-84io‘°-49H/y' 
o 


3 3 0 83  43 

where  V - 3.5  x 10  m per  kt  ’ for  events  in  alluvium,  V = 1.1  x 10  m per 
o o r 

1^0.83  fQr  eventg  jn  Coral  sand  with  a normal  device  and  = 7 x 10  m per 

kt° *83  f.Qr  tanj<  shots. 


Crater  aspect  ratios  were  found  to  be  independent  of  scaled  yield  and 
2 

the  shape  factor  V/ttR  d is  only  weakly  dependent  on  scaled  yield  (to  the  0.041 
power).  By  comparison  high  yield  Pacific  nuclear  events  exhibit  a definite  yield 
dependence  of  aspect  ratio.  Deviations  of  individual  Pacific  crater  volumes  and 
aspect  ratios  from  the  overall  trend  of  the  Pacific  data  with  yield  show  a strong 
correlation  indicative  of  either  cause  and  effect  or  mutual  cause. 

1.3  RECOMMENDATIONS 

The  centrifuge  tests  which  have  been  performed  so  far  indicate  the 
importance  of  attention  to  detail  in  test  bed  preparation.  If  care  is  not  used, 
scatter  in  the  laboratory  will  be  as  large  as  in  field  events,  and  important  features 
of  the  data  may  be  obscured.  The  tests  conducted  by  Schmidt  have  shown  both 
that  careful  preparation  can  lead  to  excellent  reproducibility  and  that  seemingly 
minor  variations  in  test  bed  assembly  can  lead  to  significant  changes  in  the 
results. 

We  suggest  that  the  next  series  of  tests  be  conducted  in  saturated  coral 
sand  in  an  effort  to  reproduce  the  results  of  the  PACE  series  of  relatively  small 
(1000  lbs.  of  TNT)  high  explosive  events.  If  the  simulation  of  these  events  proves 
successful,  an  attempt  to  reproduce  the  CACTUS  results  on  a small  scale  might 
be  undertaken  next.  This  would  require  the  development  of  a high  explosive 
equivalence  for  this  event. 

The  centrifuge  facilities  currently  available  are,  too  small  to  achieve  a 
direct  simulation  of  a megaton  nuclear  surface  burst.  In  order  to  understand  the 
large  aspect  ratios  observed  on  the  events  at  the  PPG  with  yields  above  40  kT, 
it  will  be  necessary  to  hypothesize  a physical  mechanism  to  be  tested  on  the 
centrifuge  with  a test  bed  material  modified  to  enhance  the  proposed  effect  at 
smaller  scale.  The  possibility  of  liquifaction-enhanced  slumping  is  one  example 
of  a mechanism  which  could  be  explored  through  selected  variation  in  test  bed 
properties. 

Once  the  general  trend  of  the  high  yield  data  is  understood,  one  could 
proceed  to  a detailed  study  of  each  source,  emplacement  and  geologic  media  in 
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an  attempt  to  differentiate  the  scatter  from  the  general  trend  for  individual 
events.  We  suggest  that  a return  to  geologic  media  of  strategic  interest  be 
given  higher  priority  at  this  stage  of  program  evolution. 
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2.  GEOTECHNICAL  TESTING  BY  CENTRIFUGE 


2.1  BACKGROUND 

The  use  of  a centrifuge  to  test  the  response  of  model  soil  and  rock 
structures  was  initiated  at  about  the  same  time  in  the  early  1930's  in  both  the 
United  States  by  P.  B.  Bucky, and  in  the  Soviet  Union  by  G.  I.  Pokrovsky.^  ^ 
Bucky  performed  tests  on  the  stability  of  mine  openings  by  using  very  small 
models  in  a centrifuge  having  an  arm  radius  of  about  0.2  meters.  The  tests 
were  conducted  at  centrifugal  accelerations  up  to  2,000  g's.  Bucky  postulated 
that  his  techniques  could  be  extended  to  accelerations  of  up  to  15,000  g's  on  an 
apparatus  with  an  arm  radius  of  3 meters.  Apparently,  the  technique  did  not 
meet  with  wide  acceptance  in  geotechnical  work  in  the  United  States,  and  only 
a few  subsequent  experiments  have  been  carried  out. 

In  contrast,  the  work  begun  by  Pokrovsky  and  his  colleague  I.  S.  Fyodorov 

in  the  USSR  led  to  widespread  use  of  centrifuges  for  soil  engineering  purposes 
(21)* 

in  the  Soviet  Union.  The  USSR  effort  has  also  considered  explosive  phenom- 
ena. They  have  conducted  simulations  of  explosive  loading  of  underground 
structures  and  simulations  of  HE  cratering  for  both  surface  bursts  and  shallow 
depths  of  burial/10’^ 

In  the  early  1960's,  H.  Ramberg,  a Swedish  geologist  and  a number  of 

his  students  began  an  extensive  series  of  tests  of  model  geological  structures,  in 

(22  23) 

studies  of  the  stability  of  various  tectonic  processes.  ’ In  the  middle  1960's, 
A.  N.  Schofield  in  England  investigated  the  stability  of  clay  slopes  under  rapid 
drawdown  conditions.  ' He  constructed  a large  centrifuge  having  an  arm 
radius  of  1.5  meters,  and  capable  of  carrying  a 750  kg  package  to  an  acceleration 
of  130  g's.  This  machine  was  constructed  at  the  University  of  Manchester  Institute 

of  Science  and  Technology  (UMIST).  An  even  larger  centrifuge  was  constructed 

(27) 

at  the  University  of  Manchester  under  the  supervision  of  P.  W.  Rowe. 
Experience  has  now  been  accumulated  from  a large  number  of  model  tests  of 

soil  structures  carried  out  on  these  centrifuges.  At  Cambridge  University,  where 

(28) 

Schofield  began  his  centrifuge  studies,  work  by  K.  H.  Roscoe  and  C.  P.  Wroth, 

♦Although  the  publication  date  of  this  book  is  1969,  its  English  translation  has 
only  recently  become  available  through  the  efforts  of  A.  N.  Schofield. 
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continued  on  the  centrifuge  facility  rented  from  a British  aerospace  company 
while  design  and  construction  proceeded  on  a much  larger  facility  (5  meters 
radius,  acceleration  to  300  g's).  This  centrifuge  has  been  completed  and  preliminary 
studies  are  now  under  way. 

Geotechnical  work  has  also  been  done  on  centrifuges  in  Japan,  Denmark, 
(29  33) 

France  and  South  Africa.  A literature  review  of  centrifuge  testing  in  soils 

is  given  in  Reference  34. 

In  1975,  a symposium  on  the  use  of  centrifuges  in  modeling  geotechnical 
phenomena  was  organized  at  the  California  Institute  of  Technology  by  R.  F.  Scott. 
Representatives  from  the  U.  K..  Sweden,  Canada,  and  many  U.  S.  universities 
gathered  and  presented  their  research  using  this  technique.  Representatives  of 
the  Ames  Aeronautical  Laboratory  proposed  refurbishing  their  centrifuge  for  this 
purpose.  This  would  provide  a facility  to  accelerate  payloads  up  to  approximately 
five  metric  tons  to  as  much  as  200  g's.  The  arm  radius  is  17  m and  the  workspace 

3 

would  be  10.6  m . 

Sandia  Laboratories  is  currently  operating  two  centrifuges,  with  a 7.6 

and  10.7  meter  radius  arm.  The  arm  has  a design  strength  of  7.3  x 10^g-kg  to 

200  g's.  Larger  g's  are  possible.  The  working  section  is  enclosed,  and  provides 

3 

a space  of  about  7.1  m . 

In  the  USSR,  the  cratering  due  to  buried  charges  has  been  investigated 

(21) 

for  sands  and  clays  having  various  moisture  content.  In  particular,  an 

experiment  was  conducted  at  the  VODGEO  Institute  which  utilized  a plastic  clay 
with  the  moisture  content  of  18%.  A g enhancement  of  100  was  used.  The 
Russians  claim  agreement  between  the  model  crater  and  the  prototype  crater 
volume  to  5.8%.  Tests  have  also  been  conducted  in  clayey  loam.  By  emplacing 
thin  layers  of  sand  in  these  media,  the  Soviet  researchers  have  been  able  to 
measure  the  deformation  of  the  medium  near  the  crater. 

2.2  THEORY 

In  this  section  it  will  be  shown  that  small  model  HE  tests  in  a centrifuge 
simulate  large-scale  HE  events  at  normal  gravity.  Gravity  plays  a more  important 
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role  in  the  large-scale  event,  and  the  enhanced  g's  in  the  small  experiment 
simulate  this  effect.  If  the  experiment  is  properly  designed,  many  important 
material  properties  scale  exactly,  provided  that  the  model  is  made  of  the  same 
material  as  the  prototype.  Thus,  it  is  not  necessary  to  seek  or  develop  new 
materials  to  use  in  the  model. 

It  is  shown  that  if  the  dimensions  of  the  model  are  reduced  by  the  same 
factor  as  the  g enhancement,  then  the  usual  "hydrodynamic"  scaling  laws  are 
obeyed.  This  scaling  leaves  the  units  of  velocity,  specific  internal  energy,  density, 
and  pressure  invariant.  As  a result,  material  moduli  and  their  dependence  on 
specific  internal  energy  (or  temperature)  are  preserved.  Plastic  effects  of  a very 
general  nature  are  preserved,  including  those  described  by  the  most  complex  cap 
models.  As  is  usual  in  hydrodynamic  scaling,  rate  effects  or  effects  having  an 
inherent  characteristic  length,  do  not  scale.  However,  failure  or  fracture  processes 
which  are  described  in  terms  of  strains  and/or  stresses,  will  scale  correctly. 

The  equation-of-state  of  high  explosives  and  air  scale  properly.  Thus, 
the  same  high  explosives  could  be  used  in  the  model  as  in  the  prototype.  If 
difficulties  are  encountered  with  attaining  detonation  in  small  models  of  the  HE 
used  in  the  prototype,  then  an  alternate  HE  having  a similar  Chapman-Jouget 
pressure  and  thermodynamically  similar  explosive  products  can  be  used.  Since 
the  EOS  of  air  scales  properly,  the  effects  of  air  blast  and  overpressure  will  be 
correctly  modeled. 

2.2.1  Hydrodynamic  Scaling 

It  will  be  assumed  that  readers  of  this  report  are  familiar  with  the  usual 
"hydro"  type  of  scaling.  In  "hydro"  scaling  the  physical  parameters  of  density, 
pQ,  modulus,  and  a characteristic  length,  D,  are  prescribed.  We  adopt  the 
point  of  view  that  Pq  and  W are  constant,  but  D is  permitted  to  vary.  For 
example,  it  can  be  the  characterisic  dimension  of  one  of  a set  of  geometrically 
similar  configurations.  As  shown  in  Table  1,  a characteristic  velocity  is  prescribed 
by  the  sound  speed,  c = p /p  so  that  a unit  of  time  = D/c  can  be  constructed. 
This  implies  that  as  the  unit  of  length  is  increased,  the  unit  of  time  is  increased 

♦The  equation-of-state  can  be  a non-linear  function  of  density  and  specific  internal 
energy  in  which  the  parameters  PQ  and  m0  appear. 
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Table  1 

Comparison  of  Sealings 


"HYDRO" 

P,  U,D 

FROUDE 

CENTRIFUGE 

QUANTITY 

DIMENSIONS 

c2=M/P 

P.g.  D 

2 

gD=U  -const 

Length 

L 

D 

D 

D 

Mass 

M 

pD3 

pD3 

pD3 

Time 

T 

D/c 

/D/g 

D/U 

Velocity 

L/T 

c 

/Id 

U 

Pressure 

M/LT2 

PC2 

PgD 

pu2 

Specific 

Internal 

l2/t2 

2 

C 

gD 

u2 

Energy 

Energy 

ml2/t2 

2„3 
pc  D 

PgD4 

pU2D3 

Kinematic 

Viscosity 

l2/t 

cD 

DSgD 

UD 
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by  the  same  factor.  As  Table  1 shows,  the  units  of  velocity,  pressure,  and 

specific  internal  energy  remain  invariant  under  this  type  of  size  sealing.  The 

constitutive  relations  remain  invariant  also  if  they  are  couched  only  in  terms  of 

strain,  stress,  and  specific  internal  energy.  Nonlinear  or  hysteretic  models  still 

scale  appropriately  as  long  as  they  are  couched  in  terms  of  the  above  variables. 

(35) 

This  includes  complex  void  crushing  models,  such  as  the  P-a  model  of  Herrmann. 

It  also  includes  plastie  models,  such  as  the  von  Mises,  Mohr-Coulomb,  or  Prager- 

Drucker,  and  plastic-dilatant  models  as  prescribed  by  the  cap.  It  is  seen  that 

3 1/3 

the  energy  behaves  as  D , which  implies  that  D is  proportional  to  W , which 

1/3 

is  the  genesis  of  the  widely  used  term  "W  scaling". 

The  unit  of  kinematic  viscosity  is  equal  to  cD.  This  means  that  when 
the  differential  equations  describing  the  motion  are  rendered  dimensionless  by 
this  set  of  units,  the  dimensionless  kinematic  viscosity  will  be  larger  for  smaller 
scale  models.  Thus,  the  effeets  of  viscosity  are  enhanced  in  smaller  events. 
Viscosity  is  just  one  manifestation  of  a strain  rate  effect,  and  one  can  deduce 
that  in  general  strain  rate  effects  are  enhanced  at  smaller  sizes. 

2.2.2  Froude  (Gravity)  Sealing 

If  one  considers  the  effects  of  gravity,  he  has  at  his  disposal  the  physical 
units  of  (>,  g,  and  D.  The  units  generated  from  these  quantities  are  shown  under 
the  heading  Froude  in  Table  1.  As  the  table  shows,  the  unit  of  pressure  is  given 
by  gD.  The  dependence  of  the  unit  of  pressure  on  D destroys  dynamic  similitude. 
That  is,  if  one  considers  a smaller  model,  the  dimensionless  moduli  which  appear 
in  the  constitutive  relations  will  be  larger  than  in  the  large  prototype.  This 
implies  that  in  a constant  gravity  field,  a large-scale  experiment  will  behave  as 
though  the  moduli  were  weaker  than  in  a small-scale  one.  Similarly,  the  unit 
of  specific  internal  energy  is  gD.  The  internal  energy  released  per  gram  is  an 
important  parameter  for  high  explosives,  and  the  dependence  of  this  unit  on  D 
provides  a serious  perturbation  to  scaling  high  explosive  phenomena.  It  is  seen 

4 

that  the  energy  in  this  sealing  depends  on  D , so  that  when  this  type  of  scaling 
is  appropriate,  lengths  vary  as  and  times  as 

2.2.3  Centrifuge  Scaling 

The  key  point  of  this  section  is  that  the  enhanced  effects  of  gravity  in 
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larger  events  can  be  taken  into  account,  and  dynamic  similitude  preserved,  if 

one  scales  the  unit  of  length  down  by  the  same  factor  that  the  unit  of  gravity 

is  scaled  up.  This  will  be  termed  centrifuge  scaling.  In  this  scaling,  the  product 

2 

gD  is  a constant  which  we  term  U . U has  the  dimensions  of  velocity,  and  it 
is  seen  from  the  last  column  in  Table  1 that  it  plays  the  same  role  in  determining 
units  as  c does  in  hydro  scaling.  The  basic  assumption  in  centrifuge  scaling  is 
that  rate  effects  are  not  important. 

As  a specific  example,  suppose  one  wished  to  simulate  a 500  metric  ton 
HE  sphere  detonated  half-buried  in  a soil  medium.  If  the  experiment  were 
modeled  in  a centrifuge  at  500  g's,  each  linear  dimension  in  the  model  would  be 
reduced  by  a factor  of  500.  The  HE  sphere  radius,  assuming  a loading  density 
of  1.5  g/cm  , would  be  4.30  m in  the  prototype  and  0.86  cm  in  the  model.  The 
experiment  would  start  with  the  small  sphere  half-buried  in  the  same  soil  as  the 
prototype,  and  with  air  above  it.  When  the  model  charge  is  detonated  at  its 
center,  the  explosive  detonation  wave  propagates  at  the  same  speed  as  it  would 
in  the  prototype,  and  generates  the  same  pressure.  Since  the  model  is  500  times 
smaller  than  the  prototype,  the  wave  reaches  homologous*  points  in  the  model 
at  times  l/500th  of  those  in  the  prototype.  In  particular,  it  reaches  the  sphere's 
surface  in  l/500th  of  the  time.  When  the  shock  propagates  into  the  ground  and 
air,  it  generates  shocks  having  the  same  pressures  as  in  the  prototype,  and  having 
the  same  velocities.  This  is  true  even  if  the  model  is  soil  having  extreme 
non-linearity  and  hysterisis.  The  shock  propagating  into  the  air  is  similarly 
modeled. 


Again,  the  shocks  attain  homologous  points  in  1 /500th  of  the  time  in 
the  model  because  of  its  smaller  size.  The  pressure  and  velocity  profiles  behind 
the  shocks  should  be  scale  models  of  each  other  at  corresponding  (scaled)  times. 
The  gravitational  effects  are  the  same  in  the  model  and  prototype  because  the 
500-fold  gravity  in  the  model  acts  for  1 /500th  of  the  time.  Thus,  the  ejecta 
trajectories  in  the  model  should  appear  exactly  as  those  in  the  prototype  if  slow 
motion  photography  of  the  model  slowed  them  down  a factor  of  500.  Other 
ground  motioas,  occasioned  by  both  the  direct  shock  and  by  the  air  overpressure, 
should  also  be  scale  models  of  each  other  in  both  time  and  space.  When  the 
event  is 

♦Homologous  points  are  corresponding  points  in  the  model  and  prototype. 
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over,  the  crater  and  ejecta*  in  the  model  should  be  scale  models  of  the  real 
event. 

2.3  MINIMIZING  THE  CORIOLIS  EFFECT 
. • 

The  Coriolis  acceleration  ac  in  a rotating  system  is  given  by 

_->•  -»■  . . 

ac  = -2d)  x v (5) 

—y 

where  u is  the  angular  velocity  of  the  rotating  system,  and  v is  the  velocity 
with  respect  to  the  rotating  system.  The  aeeeleration  is  produced  with  respect 
to  the  rotating  coordinate  system  by  the  tendency  of  particles  to  move  in  a 
straight  line  in  fixed  space.  Because  the  measurements  in  a centrifuge  will  be 
taken  with  respeet  to  the  rotating  system,  it  will  provide  perturbations  to  the 
experiment. 

The  magnitude  of  the  Coriolis  acceleration  is 

ac  = 2(ov  sin  y (4) 


where  Y is  the  angle  between  the  angular  acceleration  and  velocity  vectors. 
The  Coriolis  aeeeleration  is  therefore  largest  when  the  velocity  is  normal  to 
the  centrifuge's  axis  of  revolution  (Y  = 90°).  It  is  zero  when  the  velocity  is 
parallel  to  the  axis  (y  = 0°).  Thus,  there  are  two  directions  in  the  environment 
in  which  the  Coriolis  aeeeleration  is  zero. 


A measure  of  the  perturbation  is  given  by  the  relative  magnitude  of  the 
Coriolis  acceleration  and  that  of  the  centrifuge,  which  we  call  g.  Therefore,  the 
perturbation  is  given  by 


ac/9 


2(i)V  sin  y 
W2r 


2v  sin  y 
/gr 


(6) 


We  can  see  from  this  equation  that  ap/g  is  smallest  when  g and  r are 

♦This  is  true  if  Coriolis  accelerations  can  be  neglected  compared  to  the  centrifugal 
acceleration.  To  minimize  Coriolis  accelerations,  one  should  use  as  long  a 
centrifuge  arm,  and  work  at  the  highest  g's  possible.  This  is  discussed  in  the 
next  section. 
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largest.  (We  recall  from  the  discussion  on  scaling  that  the  velocity  is  invariant 
at  homologous  points  in  the  model  and  prototype.)  Therefore,  it  is  desirable  to 
work  at  the  largest  value  of  g and  r possible  to  minimize  effects  due  to  the 
Coriolis  acceleration.  If  we  work  at  a given  g level,  the  Coriolis  effect  is 
minimized  by  using  the  largest  radius  arm  possible. 

Lack  of  symmetry  in  the  final  crater  shape  provides  an  indication  when 
Coriolis  effects  have  been  important.  No  such  effect  was  observed  in  the  Boeing 
tests  discussed  in  Section  4. 
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3.  HIGH  EXPLOSIVE  CHARGE  EQUIVALENCE  FOR  JOHNIE  BOY 


A major  portion  of  our  program  was  devoted  to  determination  of  the 

spherical  PETN  charge  which  best  simulates  the  JOHNIE  BOY  event.  The  JOHNIE 

BOY  source  had  a yield^^  of  0.5  + 0.2  kT  and  a depth  of  burial  of  0.585  m. 

This  event  was  conducted  at  Area  18  of  the  Nevada  Test  Site  (NTS).  The  resulting 

(36) 

apparent  and  true  crater  profiles  are  shown'  ' in  Figure  1.  These  profiles  are 

in  good  agreement  with  the  apparent  crater  dimensions  of  18.6  m,  9.33  m and 

3 (12) 

4,113  m for  radius,  depth  and  volume  given  by  Vortman. 

Several  criteria  led  to  the  choice  of  JOHNIE  BOY  for  detailed  study. 
First,  we  were  looking  for  an  event  at  a dry,  homogeneous  site.  Second,  the  yield 
of  the  event  and  size  of  the  final  crater  should  be  on  a scale  that  could  be 
simulated  at  the  Boeing  facility.  Third,  of  the  various  nuclear  cratering  events 
at  NTS,  JOHNIE  BOY  had  been  the  most  extensively  studied. 

The  main  uncertainties  in  using  JOHNIE  BOY  as  a standard  to  validate 

centrifuge  scaling  should  be  recognized  at  the  outset  of  this  discussion.  Some  of 

these  points  will  be  brought  out  in  more  detail  in  the  following  sections.  First, 

the  details  of  source  emplacement  were  not  clearly  documented  at  the  time  of 

the  event.  Sensitivity  studies  have  not  been  performed  to  explore  what  variations 

in  energy  coupling  are  plausible.  The  alluvium  properties  of  Area  18  are  poorly 

known.  Data  from  other  sites  were  used  to  construct  the  material  response 

model.  While  PETN  may  be  considered  as  well  characterized  as  any  explosive, 

its  detonation  properties  are  quite  sensitive  to  small  amounts  of  air  voids.  The 

precise  packing  density  of  the  charges  used  in  the  Boeing  tests  was  not  known 

at  the  time  the  calculations  discussed  here  were  performed,  but  calorimeter 
(15) 

tests  on  similar  charges  suggested  that  a variation  of  ± 20%  in  detonation 
energy  might  be  expected.  Although  this  list  of  uncertainties  may  appear  formi- 
dable, we  felt  that  a simulation  which  reproduced  apparent  crater  volume  to 
within  30%  and  linear  dimensions  to  10%  was  an  ambitious,  but  not  impossible 
goal.  Results  were  22%  in  volume,  16%  in  depth  and  virtual  agreement  in  radius. 

In  Section  3.1,  we  present  the  criteria  used  to  establish  charge  equiv- 
alence. Material  models  are  summarized  in  Section  3.2.  Results  of  the  final 
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calculation  are  given  in  Section  3.3.  Sensitivity  studies  are  summarized  in  Section 
3.4.  The  results  of  the  centrifuge  tests  are  discussed  in  Section  3.5. 

3.1  EQUIVALENCE  CRITERIA 

The  calculations  described  in  Reference  36  were  chosen  as  the  standard 
for  determination  of  charge  equivalence.  Thus,  we  introduce  the  additional 
uncertainty  of  the  effect  of  different  codes  used  by  different  groups  at  different 
times.  Comparisons  of  codes  in  the  past  using  identical  starting  conditions  have 
led  to  variations  in  the  final  results  as  large  as  any  of  the  uncertainties  mentioned 
in  the  previous  section. 

The  energy  partition  in  the  ground  beyond  a range  of  3.6  meters  is  shown 

in  Figure  2.  One  can  see  that  the  energy  coupling  is  essentially  complete  after 

several  miliseconds.  It  is  on  this  time  scale  that  we  presume  the  details  of  the 

source  become  unimportant  and  the  nuclear  to  high  explosive  equivalence  can  be 

(37) 

established.  As  in  previous  studies  for  buried  cratering  sources,  attention  is 
initially  focused  on  the  kinetic  energy  in  the  flow. 

(38 ) 

Maxwell  noted  from  a study  of  several  surface-burst  calculations  that 
the  cratering  flow  had  features  consistent  with  steady-state  flow  of  an  incom- 
pressible fluid  in  the  region  behind  the  shock  wave  after  the  stresses  had  relaxed 
to  low  values.  In  order  to  achieve  a quantitative  measure  of  this  observation, 
concentric  surfaces  S ^(0)  (i  = 1,2,...)  were  prescribed  at  t = 0 as  illustrated  in 
Figure  3a.  Each  surface  was  defined  by  Lagrange  grid  points  in  the  computational 
mesh  that  moved  with  the  flow  as  illustrated  in  Figure  3b.  The  volume  V.(t) 
enclosed  by  the  surface  S.(t)  and  the  original  ground  plane  and  its  time  rate  of 
change  V.(t)  were  observed  to  obey  the  relationship 

v.(t)  = e/v^  (t)  (7) 

where  i and  P were  independent  of  time  and  index  i after  the  shock  had  passed 
and  the  stresses  had  relaxed.  If  one  defines  the  mean  radius 

H = (3V/2tt}1/3,  (8) 
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DEPTH  (m) 


RADIUS  (m) 


Figure  1.  Two  sections  of  JOHNIE  BOY  apparent 
crater  (solid  lines)  and  true  crater 
(dashed  lines). 
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Partition  of  energy  versus  time  in  the  ground  beyond  the 
range  of  3. 6 meters. 


Figure  3b. 


& 


Reference  surface  and  its  associated 
volume  at  t=0  and  t>0. 


i 
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and  mean  radial  velocity 

R = V/2t,R2,  (9) 

the  equivalent  observation  is 

R = a/RZ  . (10) 

One  can  think  of  a as  a measure  of  the  kinetic  energy  of  the  flow  and  Z as  a 
measure  of  the  geometry  of  the  streamlines.  For  purely  radial  flow,  Z = 2.0. 
These  general  observations  form  the  basis  for  the  hypothesis  that  a complex,  high 
energy  source  might  be  simulated  by  variation  in  the  yield  and  depth  of  burial 
of  a spherical  high  explosive  charge.  Although  the  parameters  are  cross-coupled, 
it  was  expected  that  a would  depend  strongly  on  charge  yield  and  Z would 
depend  most  strongly  on  depth  of  burial. 

(18) 

Application  of  this  analysis  to  the  JOHN1E  BOY  calculation  led  to 
the  results  shown  in  Figures  4 and  5.  One  can  see  that  the  results  are  quite 
consistent  with  Eq.  (10)  for  the  values  Z = 2.71  and  a = 2.44  x 104  in  mks 
units.  Some  deviation  is  apparent  at  times  beyond  100  msec  when  gravity  and 
strength  effects  cause  a breakdown  of  the  steady-state  assumption.  The  most 
severe  test  of  Eq.  (10)  is  illustrated  in  Figure  6,  where  the  space-time  motion 
of  the  surface  which  started  at  3.6  m is  compared  with  the  model. 

3.2  MATERIAL  PROPERTIES 

3.2.1  Alluvium 

The  material  properties  at  the  JOHNIE  BOY  site  (Area  18  at  NTS)  are 
poorly  known.  Data  from  other  sites  were  used  to  establish  the  model  for  the 
previous  calculations  of  this,  event.  We  have  used  this  same  data  base  to  choose 
parameters  in  our  model  for  the  response  of  geologic  materials.  The  initial  density, 

3 

bulk  modulus  and  shear  modulus  were  given  the  values  of  1.8  g/cm  , 7.717  kbar 
and  5.5  kbar  used  in  the  earlier  models.  The  load  and  unload  curves  which  form 
the  envelope  of  the  irreversible  compaction  region  are  shown  in  Figure  7.  An 
air-void  content  of  19%  is  presumed. 

(39) 

The  EOS  of  the  void-free  material  has  been  given  previously  by  Allen. 
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t=0.351  msec- — (J> 


600  msec 


0 5 10  15 


HORIZONTAL  RANGE  (m) 

Figure  6.  Model  velocities  (dashed  arrows)  compared 
to  computed  velocities  (solid  arrows)  at 
selected  points  and  times. 
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For  condensed  states  (p/p0>1)  and  for  any  cold  states  (e<e$) 

P = P = Ap  + Bp?  + Tpc  (11) 

where  p = ri-1  , n = p/ PQ  = vq/v 

r = a + b (12) 

c + 1 

" -JL 


For  expanded  hot  states  (p/p  <1  and  e>e') 


P = Pe  = aep  + 


bf  p _ + Ape 

e + 1 


-3(v/vq-1  )' 


e-a(v/v0-l)‘ 


(13) 


A smooth  transition  between  the  condensed  and  expanded  states  is  insured 

by  a transition  equation  for  the  intermediate  region  defined  by  e,.<c  < c ' and 

b s 

p/P0<l.  This  blended  portion  of  the  equation  of  state  has  the  form: 


P = 


(c-cs)  Pe  + (e'-e)  Pc 
e_s  - cs 


(14) 


In  these  eq-  ations,  e,  p and  P are  density,  specific  internal  energy,  and  pressure. 
The  values  jf  the  constants  used  for  JOHNIE  BOY  alluvium  are  given  in  Table 
2. 


A non-assoeiative  flow  rule  was  used  with  the  yield  surface  shown  in 
Figure  8.  At  the  early  times  of  flow  development  studied  here,  the  details  of 
the  treatment  of  the  strength  of  the  alluvium  are  unimportant. 

3.2.2  PETN 

The  Jones-Wilkins-Lee  (JWL)  equation  of  state  for  detonation  products 
of  explosives  can  be  written  in  the  form: 

-R,v  ..  -R9v 

p = A(l-  — -)e  + B(  1-  p~)e  + wpe  . (15) 
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Table  2 

Parameters  for  J0HN1E  BOY  Alluvium 


Parameter 

Value 

Units 

A 

2.34 

X 

1011 

dynes/cm 

B 

1.02 

X 

1011 

dynes/cm 

a 

0.5 

- 

b 

1.3 

- 

po 

2.22 

g/cm3 

Es 

3.5 

X 

1010 

ergs/g 

cs 

1.8 

X 

io11 

ergs/g 

£ 

0 

2.01 

X 

108 

ergs/g 

a 

5.0 

- 

6 

5.0 

- 

30 


50 

40 

30 

o o o 

CNJ  r— 

(s^eq)  tcf 

l 

31 

where  c is  the  specific  internal  energy,  p is  density  and  v = P /p  . The  values 

3 U 

of  the  various  parameters  for  PETN  (with  PQ  = 1.77  g/em  ) are  given  in  Table  3 
along  with  the  heat  of  detonation  and  detonation  velocity  D.  As  was  mentioned 
earlier  in  this  section,  the  detonation  properties  of  PETN  are  sensitive  to  peeking 
density.  This  variability  is  illustrated  in  Figure  9 from  Reference  ".  Calorimeter 
tests  indicated  a variability  of  up  to  20%  in  heat  of  detonation  for  barges  similar 
to  those  used  in  the  Boeing  tests. 

3.3  RESULTS  OF  FINAL  SIMULATION  PREDICTION 

A number  of  eharge  sizes  and  depths  of  burial  were  investigated  before 
our  final  ehoiee  was  made.  The  results  of  these  preliminary  calculations  are  briefly 
summarized  in  Seetion  3.4.  Our  last  calculation  incorporated  a eharge  with  an 
initial  radius  of  1.88  m and  a depth  of  burial  of  1.20  m.  This  eharge  has  a total 
mass  of  4.93  x 104  kg  and  a total  yield  of  2.81  x 1011  J.  Thus,  the  high  explosive 
to  nuclear  energy  equivalence  faetor  represented  by  this  eharge  is  0.134.  ' 

The  initial  grid  configuration  is  shown  in  Figure  10.  An  expanded  plot  of 
the  zoning  in  the  high  explosive  sphere  and  its  immediate  surroundings  is  given  in 
Figure  11.  The  first  0.2  msee  of  the  calculation,  which  is  the  time  for  the 
spherical  HE  bum,  was  performed  with  a one-dimensional  code  so  that  the  develop- 
ment of  the  detonation  wave  eould  be  provided  adequate  resolution.  Profiles  of 
pressure  and  partiele  velocity  at  the  time  of  overlay  into  the  two-dimensional  grid 
are  shown  in  Figures  12  and  13.  The  final  grid  configuration  in  the  ground  and 
lower  portion  of  the  high  explosive  is  shown  in  Figure  14  while  the  interfaces  in 
the  upper  region  of  interaction  with  the  air  are  illustrated  in  Figure  13. 

As  ean  be  seen  in  Figure  11,  a hemi-spherieal  surface  in  the  ground  at  a 
radius  of  3.6  m was  established  in  the  initial  grid  generation  and  maintained  for 
editing  purposes  and  comparison  with  the  equivalence  criteria  outlined  in  Seetion 
3.1.  The  kinetic  energy  in  ground  material  beyond  this  surfaee  is  compared  with 
the  JOHNIE  BOY  results  in  Figure  16.  The  motion  of  the  reference  surfaee, 
averaged  as  diseussed  in  Seetion  3.1,  is  compared  with  the  JOHNIE  BOY  results  in 
Figure  17.  We  consider  this  comparison  to  be  very  favorable.  Finally,  the  detailed 
velocity  fields  of  the  reference  surfaee  are  compared  in  Figure  18.  Again  the 
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Table  3 


Parameter 


EOS  Parameters  for  PETN 
Value 


Units 


po 

1.77 

g/cm3 

R1 

4.4 

- 

R0 

1.2 

- 

c, 

A 

6.17 

X 

1012 

dynes/cm 

B 

1.69 

X 

io11 

dynes/em 

U) 

0.25 

t— * 

o o 

CO  *— * 

o 

- 

Go 

D 

5.7 

8.3 

X 

X 

ergs/g 

m/s 

36 


33 


30 


27 


i 24 


21 


18 


15 


12 
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PRESSURE  AT  CHAPMAN- UOUGET  PLANE  (DETONATION  PRESSURE )-atm  x 10 


RADIUS  (m) 

Figure  12.  Particle  velocity  distribution  in  PETN  at  termination  of  burn. 


Figu 


RADIUS  (m) 


Figure  14.  Final  grid  configuration  near  bottom  half 
of  the  source. 
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HEIGHT  (m) 


0 10.1 


RADIUS  (m) 


Figure  15.  Final  grid  configuration  and  material 
interfaces  above  the  source. 
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Figure  16.  Kinetic  energy  coupling  for  final  prediction. 


R (m/s) 


DEPTH  (m) 


comparison  appears  quite  good. 


The  final  step  in  defining  the  simulation  event  on  the  Boeing  facility  is 
to  choose  the  gravitational  acceleration  to  be  used  with  one  of  the  standard 
PETN  charges  that  were  available.  The  1.2  gm  charge  was  chosen  which  led  to  a 
value  of  g given  by: 


g = 


7 

4.93  x 10  gm 
1.2  gm 


1/3 


345, 


and  the  depth  of  burial  ^should  be  120/345  = 0.348  cm. 


3.4  SENSITIVITY  STUDIES 

Several  additional  calculations  were  run  for  a few  miliscconds  as  we 
iterated  on  both  high  explosive  yield  and  depth  of  burial  to  determine  the  best 
simulation  charge.  In  the  first  series,  the  depth  of  burial  of  the  high  explosive 
was  fixed  at  the  center  of  the  JOHNIE  BOY  source  emplacement.  Calculations 
were  made  with  initial  charge  radii  of  1.75;  2.2  and  2.57  m.  The  kinetic  energy 
coupling  beyond  the  3.6  m reference  surface  is  shown  in  Figure  19.  The  results 
indicate  that  the  coupled  energy  increases  more  rapidly  than  a simple  linear 
function  of  charge  yield. 

When  wc  looked  in  more  detail  at  the  motion  of  the  reference  surface, 
it  became  clear  that  an  increased  burial  depth  for  the  high  explosive  would 
improve  the  simulation.  The  results  of  the  two  calculations  at  a depth  of  1.2  m 
are  shown  in  Figure  20.  While  these  results  suggest  that  the  high  explosive  yield 
can  be  fine  tuned  to  approximately  ±,10%,  it  is  not  as  clear  how  well  the  "best" 
depth  of  burial  can  be  determined.  This  choice  hinges  more  on  flow  field 
comparisons  such  as  those  shown  in  Figure  18.  We  doubt  that  an  increase  in 
charge  burial  will  improve  this  comparison  to  any  marked  degree. 

3.5  COMPARISON  WITH  SIMULATION  TEST  RESULTS 

The  dimensions  of  the  predicted  simulation  of  JOHNIE  BOY  are  compared 
with  the  Boeing  test  results  in  Tabic  4.  The  test  JB1  was  conducted  at  the 
suggested  depth  of  burial  of  0.348  cm  (0.64  charge  radii  for  the  1.2  gm  charge) 
as  accurately  as  possible.  As  a study  of  depth  of  burial  sensitivity,  the  DOB 
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Figure  19.  Effect  of  PETN  charge  size  at  0.53  m DOB. 


Table  4 

Comparison  of  Simulation  Crater  Dimensions 


JOHNIE 

BOY* 

JB1 

JB2 

Radius  (cm) 

5.39 

5.37 

6.08 

Depth  (cm) 

2.70 

2.26 

2.76 

Apparent  ^ 

Volume  (citi 

100. 

77.7 

116. 

Aspect  Ratio 

2.0 

2.38 

2.2 

♦Scaled  down  by  345 

in  linear  dimensions. 
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n ww 


was  increased  to  1.47  charge  radii  for  JB2.  The  effects  of  DOB  variation  are 
discussed  in  more  detail  in  Section  4.  We  consider  the  simulation  quite  satisfactory 
considering  the  uncertainties  we  have  already  discussed. 

We  would  like  to  emphasize  once  more  the  uncertainties  which  are  likely 

to  frustrate  further  refinement  of  this  simulation.  Measurements  of  density  at 

3 

the  JOHNIE  BOY  site  have  varied  from  a low  of  1.39  g/cm  to  a high  of  1.8 

3 

g/cm  and  moisture  content  has  shown  a variation  from  7.1  to  11.5%.  It  has 
(IB) 

been  suggested  ' that  the  details  of  source  emplacement  can  change  the  charge 
equivalence  by  30%.  The  yield  of  the  PETN  charges,  which  wc  have  suggested 
as  a third  source  of  uncertainty,  could  be  better  determined  by  further  testing. 
Such  tests  should  be  conducted  if  only  to  establish  the  true  TNT  equivalence  of 
these  particular  charges.  In  addition,  the  observed  variation  could  be  caused  by 
a weak  rate  effect  causing  a deviation  from  true  similitude  in  these  small-scale 
tests. 
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4.  ANALYSIS  OF  CENTRIFUGE  CRATERING  DATA 


The  purpose  of  this  section  is  to  analyze  the  centrifuge  crater  data  obtained 
by  Boeing  in  conjunction  with  preexisting  data  obtained  at  1 g.  We  describe  several 
systematic  relations  between  dimensions  of  explosive  craters  produced  by  near- 
surface explosions  and  parameters  of  the  explosion,  such  as  yield  and  height  of 
burst.  This  comparison  demonstrates  the  power  of  the  centrifuge  approach,  as  well 
as  providing  potentially  useful  relations  in  its  own  right.  Some  of  the  systematic 
relations  are  expected  from  theoretical  considerations,  but  others  are  purely  empir- 
ical. It  is  hoped  that  the  delineation  of  the  empirical  relations  will  stimulate 
theoretical  analysis  of  the  cratering  process  in  order  that  their  bases  may  be 
ascerta  ined. 

4.1  BACKGROUND 

(41) 

Earlier  work  by  Cooper  has  related  the  effect  of  height  of  burst  on 

cratering  efficiency  (crater  volume  divided  by  explosive  yield),  and  has  achieved 

considerable  success  in  synthesizing  cratering  data  from  a wide  range  of  yields, 

explosive  types  and  geologic  media.  Briefly,  Cooper  suggested  that  the  apparent 

crater  volume  V produced  by  an  explosion  of  yield  Y could  be  expressed  as 
8 

V = v(G,S,H)  x Y (16) 

8 

where  the  cratering  efficiency  v(G,S,H)  is  a function  of  the  geology  (G),  the  explosive 

source  (S)  and  the  scaled  height  of  burst  (H).  The  parameter  used  to  scale  the 

1/3  1/3 

height  of  burst  and  other  linear  dimensions  was  chosen  to  be  V (i.e.,  H=h/V  ; 

8 8 

h = height  of  burst).  Cooper  factors  v(G,S,H)  into 

v(G,S,H)  = K(G)vq(Go,S)F(S,H)  (17) 

where  vo  (Gq,  S)  is  a "standard"  cratering  efficiency  for  source  S at  zero  height 
of  burst  in  a standard  geology  G , K(G)  is  the  ratio  of  the  zero  height  of  burst 
efficiency  in  geology  G to  that  in  Gq,  and  F(S,H)  is  a height  of  burst  shape  factor 
for  the  source  S. 

Recently  two  series  of  very  precisely  controlled  cratering  experiments 
have  been  conducted  which,  when  considered  together,  indicate  a fundamental 
difficulty  in  the  volume  scaling  analysis.  Results  of  these  experiments  are 
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(15) 

summarized  in  Tables  5 and  6.  The  first  set  of  data,  obtained  by  Piekutowski, 
consist  of  a large  number  of  cratering  experiments  in  Ottawa  sand  using  small 
(0.4  gm  to  4 gm)  high  explosive  spherical  charges.  Considerable  effort  was 
directed  to  producing  precisely  characterized  and  reproducible  test  beds.  As  a 
result,  the  data  obtained  were  very  reproducible.  Piekutowski  investigated  the 
effects  of  charge  type  (PETN  and  PbfN^  were  used),  height  of  burst  and  test 
bed  bulk  density.  Observed  cratering  efficiencies  were  an  order  of  magnitude 
higher  than  efficiencies  obtained  with  100  kg  to  1000  kg  explosive  charges  in 
natural  sands  such  as  alluvium/42^ 

(43) 

Schmidt  has  conducted  tests  similar  to  those  of  Piekutowski,  using 
0.4  to  4 gm  charges  in  Ottawa  sand  at  zero  height  of  burst,  but  at  high  acceleration 
in  a centrifuge.  His  data,  together  with  Piekutowski’s,  show  that  cratering 

3 

efficiency  decreases  as  the  product  g Y increases.  Specifically, 

vq  a Y5/6  //g  . (18) 

Since  yields  of  the  Pacific  cratering  events  are  typically  3 or  4 orders  of  magnitude 
larger  than  those  of  Nevada  events,  application  of  Eq.  (18)  would  indicate  standard 
cratering  efficiencies  lower  by  a factor  of  10  to  10  (3.2  to  4.6)  ; hence, 

the  data  base  was  re-examined  to  see  if  use  of  Eq.  (18)  instead  of  Eq.  (16)  could 
lead  to  an  improved  predictive  capability.  Schmidt  has  also  conducted  a height 
of  burst  study  at  high  g-levels  in  a reconstituted  dry  alluvium.  These  data  are 
discussed  in  Section  4.3  of  this  report. 


4.2  SCALING  RELATIONS 


Numerous  scaling  relations  have  been  proposed  for  reconciling  crater 
dimensions  from  tests  at  different  yields  and  in  different  materials.  A compre- 
hensive review  of  them  is  not  the  object  of  thL  eport.  The  interested  reader 

(2-9  12) 

is  referred  to  the  many  published  papers  dealing  with  the  subject.  ’ The 

(44) 

brief  discussion  here  is  based  on  a dimensional  analysis  by  Schmidt.  He  groups 
the  variables  into  four  dimensionless  groups  as  follows: 


U1 

7T2 

*3 

^4 


V W"1  p 

g W1'3  U-2  6-1/3 

,J/3  u-1/3 
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Table 

6.  Summary  of  UQRI  and  BAC  Cratering 

Oata  in 

Ottawa  : 

Sand 

Test 

Oensity 

Charge 

Charge 

TNT 

1 

H08 

Volume 

Radius 

Depth 

Accel . 

Type 

Mass 

Equivalent 

(Mg/m3) 

(Rg) 

(Rg) 

(m) 

(m3) 

(m) 

(m) 

(g ' s ) 

UORI-9 

1.80 

PbN. 

0 

1.7  x 10'3 

5.7 

x 10"4 

0 

4.36 

x 10*4 

1.02  x 

10'1 

2.29  x 

10'2 

i 

'JORI-22 

1.80 

Pb-'fi 

1.7  x 10*3 

5.7 

x 10‘4 

0 

4.19 

x 10'4 

1.01  x 

10'1 

2.35  x 

10"2 

1 

UDRI-107 

1.80 

petn 

4 x 10'4 

6 

x 10'4 

0 

2.59 

x 10'4 

8.73  x 

IQ'2 

1.97  x 

10'2 

i 

JDRI-1 17 

1.81 

PETN 

4 x 10‘4 

6 

x 10"4 

0 

2.74 

x 10'4 

8.76  x 

10‘2 

2.06  x 

10'2 

l 

uor:-120 

1.79 

PETH 

4 x 10'4 

6 

x 10‘4 

0 

2.92 

x 10'4 

8.88  x 

10-2 

2.13  x 

10*2 

1 

^cri-146 

1.81 

PETN 

4 x 10'4 

6 

x 10'4 

0 

2.84 

x 10"4 

8.89  x 

10'2 

2.02  x 

lO'2 

l 

SAC-10-0 

1.73 

PETH 

4 x 10'4 

6 

x 10‘4 

0 

1.71 

x 10‘5 

3.56  x 

10"2 

7.1  x 

ID'3 

463 

SAC-IO-X 

1.77 

Pb'<6 

1.7  x 10'3 

5.7 

x 10"4 

0 

2.71 

x 10'5 

3.87  x 

10"2 

8.9  x 

10*3 

463 

eAc-ii-o 

1.78 

PETM 

1.25  x 10'3 

1.7 

x 10’3 

0 

3.80 

x 10'5 

4.31  x 

10-2 

9.7  x 

10'3 

451 

8AC-11-X 

1.78 

PETN 

1.25  x 10'3 

1.7 

x 10'3 

0 

3.51 

x ,10’5 

4.38  x 

10'2 

9.9  x 

10*3 

451 

BAC- 12-0 

1.78 

PETH 

4 x 10"3 

5.5 

x IQ’3 

0 

1.01 

x 10‘4 

6.41  x 

10'2 

1.31  x 

10*2 

451 

BAC-12-X  * 

1.63 

PETN 

4 x 10'3 

5.5 

x 10'3 

0 

9.63 

x 10'5 

6.33  x 

10'2 

1.23  x 

10'2 

451 

BAC-13-0 

1.73 

PETN 

4 x 10'3 

5.5 

x ID’3 

0 

1.25 

x 10'4 

6.88  x 

10'2 

1.37  x 

10-2 

306 

BPC-13-X  * 

1.68 

PETH 

4 x 10'3 

5.5 

x 10'3 

0 

1.19 

x I0"4 

6.65  x 

10‘2 

1.56  x 

10-2 

306 

8AC-15-0 

1.78 

PETH 

4 x 10"4 

5.5 

x IQ*4 

0 

9.52 

x ID'5 

6.40  x 

10'2 

1.27  x 

10‘2 

10 

BAC-15-X 

1.78 

PETH 

4 x 10'4 

5.5 

x 10'4 

0 

9.48 

x IQ*5 

6.40  X 

10'2 

1.26  x 

10*2 

10 

* Two  high  g tests 

used  a finer  grain  siye 

fraction  of  Ottawa 

sand  which  did 

not  pack  as 

tightly  as  the 

coarser 

material . 
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where 


V = crater  volume 
o=  medium  density 
W = mass  of  charge* 
g = gravitational  acceleration 
U =Chapman-Jouget  particle  velocity 
6 =Chapman-Jouget  density 
h = height  of  burst. 


The  crater  volume  then  is  given  by  solution  of  f^TT-j  ,-n^  ,71^  ,tt^)  = 0-  For  almost 
all  materials,  Schmidt  considered  that  neither  p nor  <$  varies  by  an  appreciable 
amount,  so  that  we  can  solve  f^(  71 7 »TT2  *^3  ^ = For  Part>cular  case  °f 
zero  height  of  burst,  ^3  = constant  and  we  have  7T-|7T2°t  = constant 


or 


-a/  3 


(19) 


Based  on  this  analysis  we  expect  that  a plot  of  the  log  of  the  crater 

3 ri 

efficiency  (V/W)  against  the  log  of  g W will  be  a straight  line  of  slope  "<*/ j for 
tests  with  a given  explosive.  Data  for  dry  Ottawa  sand  at  lg  and  at  high  g are 
shown  in  Figure  21.  Over  a range  of  almost  nine  orders  of  magnitude  the  value 
of  a is  determined  to  be  about  1/2  when  only  data  using  PETN  are  considered. 
(The  best  fit  is  obtained  with  a =0.478.) 

C /£? 

(It  should  be  noted  that  the  a =1/2  value  indicates  that  V/Y  , rather 
than  V/Y,  should  be  used  when  comparing  craters  produced  by  quite  different 
yields.  Thus,  Eq.  (16)  should  be  modified  to  be 

Vq  = v (G,S,H)  *Yn  (20) 

where  V is  the  efficiency  in  units  of  volume/(yield)n  instead  of  Cooper's  volume/- 
a c /g 

yield  and  n=5/6.  We  shall  return  to  a discussion  of  this  Y dependence  when 
discussing  height  of  burst  effects  in  Section  4.3) 

♦When  considering  a single  type  of  explosive,  W and  Y are  interchangeable. 
However,  such  is  not  the  case  when  comparing  different  explosives. 
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The  effect  of  changing  the  type  of  explosive  is  also  readily  seen  in  Figure  21. 

Charges  of  PETN  used  actually  contain  over  0.1  gram  of  AgN^  as  an  initiator.  This 

has  little  effect  on  4 gram  and  1 gram  charges,  but  constitutes  an  appreciable  amount 

3 s 

of  the  charges  equivalent  to  0.4  grams  of  PETN.  The  high  g data  at  g Y = 10  (kg^,j, 
g's3)  were  obtained  using  the  same  0.4  gram  PETN  charges  and  are  seen  to  lie 
significantly  off  of  the  line  of  the  other  PETN  data.  In  addition,  both  high  g and  low 
g tests  were  conducted  using  PWN,^  charges  of  the  same  TNT  equivalence  as  the 
PETN  charges.  PbCN^  is  apparently  50%  more  efficient  in  cratering  than  is  PETN 
for  a given  energy  content. 


In  order  to  quantify  the  effects  of  changing  explosive  type,  we  can  look  at  the 

Ot 

relation  tt-j  = constant.  In  particular  the  ^ term  contains  information  on  the 

C-J  state  of  the  explosive.  So  if  we  plot  ir,  versus  tt?  , instead  of  just  V/Y  versus 

3 c 

g Y,  we  would  expect  data  for  all  explosives  to  lie  on  a single  line.  Using  the  explosive 

properties  in  Table  7,  we  find  that  this  is  indeed  the  case  for  Ottawa  sand,  as  illustrated 

in  Figure  22.  Hence,  by  using  the  entire  tt9  term,  we  can  bring  all  the  data  for 

3 *- 

Ottawa  sand  at  a density  of  1.78-1.80  Mg/m  onto  a single  line.  (The  effect  of  AgN^ 
initiators  is  still  seen  because  the  initiator  was  not  included  in  describing  the  C-J 
state.) 


(15  42  43) 

Cratering  data  from  alluvium,  playa  sand  and  other  soil  media  ’ ’ are 

presented  in  Figure  23  in  the  form  of  a 1T-|-7T2  pl°t*  Included  are  Schmidt's  data 
obtained  at  300  to  450  g's  and  Piekutowski's  at  1 g in  reconstituted  Kirtland  AFB 
alluvium.  With  the  exception  of  the  PACE  tests,  a series  of  1000  lb.  shots  in  saturated, 
unconsolidated  coral  sand,  the  data  can  be  fit  to  a single  straight  line  with  a standard 
deviation  of  only  36%  in  crater  volume.  The  line  shown  in  the  figure  has  the  equation 


(0.062+0.022  m3)W°’977g"^’070p-1^ 


t T° 

TNT  TNT 


j 0.023 


(21) 


3 

where  W is  the  charge  mass  in  kg,  g is  in  g's,  p is  in  Mg/m  and  U and  6 are  in 
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Figure  23.  tti-tt2  plot  for  craters  in  alluvium,  soil  and  saturated  coral  sand, 
best  fit  to  data  for  Ottawa  sand  shown  in  Figure  22. 


any  consistent  units.  These  tests  encompass  charges  from  less  than  a gram  to 

3 

40  tons  and  media  with  densities  from  1.13  to  1.70  Mg/m  . In  view  of  this  wide 
diversity,  the  narrow  scatter  of  the  data  is  remarkable. 

All  of  the  soils  that  fit  Eq.  (21)  are  at  least  moderately  cohesive.  It 
is,  therefore,  possible  that  the  difference  in  behavior  of  these  soils  from  that  of 
Ottawa  sand  is  due  to  their  cohesive  nature  which  inhibits  cratering  at  low  yields. 

e. 

It  is  perhaps  significant  that  the  cohesion  of  the  KAFB  alluvium,  about  1.5x10  Pa 
(43,  App.E),  is  approximately  equal  to  the  initial  overburden  stress  (Pgd)  at  the 

5 

bottom  of  the  crater,  1-1.4x10  Pa.  We  expect  that  at  still  larger  scaled  yields, 
the  effect  of  cohesion  would  disappear  and  crater  trends  should  parallel  those  of 
Ottawa  sand.  For  more  cohesive  materials  this  change  should  occur  at  even 
higher  yields. 

In  view  of  our  ultimate  goal  of  understanding  kiloton  and  megaton  nuclear 

craters,  we  feel  it  is  necessary  to  speculate  on  the  trend  which  might  be  followed 

by  craters  in  alluvium  at  still  higher  yields.  It  is  reasonable  to  expect  that 

cohesion  could  decrease  cratering  efficiency  since  more  work  must  be  done  to 

expand  a cohesive  material  than  a non-cohesive  one.  However,  it  is  difficult  to 

see  how  cratering  efficiency  could  ever  be  higher  in  a cohesive  material  than  in 

an  otherwise  similar  non-cohesive  material.  That  being  the  case,  we  should 

“I  /6 

expect  that  at  high  yields  the  cohesive  material  would  exhibit  the  same  Y 

dependence  of  cratering  efficiency  on  yield  as  is  seen  in  the  non-cohesive  medium. 

However,  the  existing  data  are  insufficient  to  determine  where  this  change  from 
_0  02  *1  /6 

Y to  Y should  occur.  The  presently  available  tests  extend  only  up  to 
the  intersection  with  the  Ottawa  sand  trend.  Unfortunately,  that  is  not  necessarily 
an  upper  limit.  The  PACE  tests  indicate  that  higher  efficiencies  are  possible  in 
wet  media.  A centrifuge  test  with  a 1.7  gram  PUN^  charge  at  500  g's  would 
be  most  helpful  in  this  regard  since  it  would  have  equivalent  to  about  5 x 10 
kg  of  TNT. 

If  the  above  comments  are  valid,  it  will  be  difficult  to  apply  a single 

factoring  approach  as  used  by  Cooper  (e.g.,  Eq.  (17))  over  the  entire  range  of 

— 1 /fi 

yields,  but  at  high  yields,  if  all  media  trend  as  Y , then  the  factoring  approach 
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should  apply  and  Eq.  (16)  can  be  rewritten 

V = K(G)v  (S)F(H,S)n(Y)  (22) 

where  K(GM/p,  vQ(S)  = (U61/6)/(e5/6UTNT  61/6TN?),  and  n(Y)  = Y5/6/g.  F(S,H) 
is  a height  of  burst  function  discussed  in  Section  4.3,  and  e is  the  ratio  of  the 
specific  energy  to  that  of  TNT. 

It  is  difficult  to  apply  Eq.  (21)  to  nuclear  cratering  events  since  the  C-J 

state  cannot  be  defined.  However,  noting  the  Eq.  (22)  factors,  we  can  still  plot 

the  nuclear  data  as  though  it  were  an  energy  equivalent  mass  of  TNT.  In  Figure 

23,  we  have  indicated  by  a N the  interpolated,  zero  height  of  burst,  l kt  nuclear 
(12) 

crater  of  Vortman.  Based  on  this  point,  we  see  that  vq(S)  ~ 5.0  relative  to 
the  Ottawa  sand  line  or,  as  an  inferred  upper  limit  vq(S)  ~ 7.0  relative  to  the 
low  yield  alluvium  line.  This  compares  with  a value  of  7.35  obtained  in  the 
previous  section  of  this  report  by  entirely  different  means.  We  will  return  to 
this  point  after  discussing  height  of  burst  effects. 

4.3  HEIGHT  OF  BURST  EFFECTS 

In  order  to  discuss  height  of  burst  effects  it  is  necessary  to  have  a 
rational  method  to  adjust  heights  for  tests  at  different  yields.  Since  the  dimensions 
of  volume  are  the  cube  of  the  dimensions  of  length,  Cooper  used  the  cube  root 
of  the  crater  volume  to  scale  height  of  burst.  One  result  of  such  a scheme  is 
shown  in  Figure  24.  (The  raw  data  from  nuclear  tests  Nevada,  plotted  as  X, 
have  been  adjusted  by  a factor  of  3.4  to  the  points  □ to  account  for  differences 
in  the  media.)  Cooper's  interpretation  is  that  the  tank  shots  behaved  in  a manner 
similar  to  the  "low  energy  density"  Nevada  events  whereas  the  other  shots  exhibited 
a markedly  different  height  of  burst  effect  because  they  were  "high  energy 
density"  sources.  (The  ZUNI  crater  ( a ),  which  partially  overlapped  earlier  craters, 
is  anomalous.)  The  slope  of  the  "high  energy  density"  curve  far  above  the  interface 
is  nearly  identical  to  that  of  the  low  energy  density  curve.  There  is  no  data 
base  on  which  details  of  the  curvature  of  the  high  energy  curve  can  be  evaluated. 
Hence,  it  is  difficult  to  use  this  scheme  to  correct  the  high  yield  Pacific  events 
to  a zero  height  of  burst. 

5 /6 

We  saw  in  Section  4.2  that  the  crater  volume  in  sand  varies  as  Y . 
That  being  the  case,  it  L reasonable  to  replace  by  (y5/6)l/3  _ y5  as 
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the  parameter  for  scaling  linear  dimensions.  Thus,  the  scaled  height  of  burst 


would  be  HOB/Y" 


If  we  plot  v/Y^  vs.  HOB/Y3^3,  we  should  be  able  to 


determine  the  effect  of  height  of  burst  on  crater  size  for  craters  from  a wide 

range  of  yields,  as  in  Figure  25.  Here  we  see  that  the  "high  energy  density" 

events  fall  on  the  same  curve  as  the  Nevada  events.  (The  line  is  a least  squares 

fit  to  the  eight  "high  energy  density"  points  without  the  Nevada  data.)  As  before, 

the  two  tank  shots  are  anomalous,  but  they  do  represent  the  same  slope  as  the 

rest  of  the  data.  The  scatter  of  data  about  the  curve  in  Figure  25  is  greater 

than  in  Figure  24,  but  the  fact  that  all  normal  nuclear  sources  (i.e.,  not  tank 

shots)  fit  a single  curve  is  a decided  advantage.  The  deviations  are,  to  some 

extent,  systematic  as  will  be  discussed  below.  The  greatest  advantage  of  using 

Figure  25  instead  of  Figure  24  is  that  the  intercept  of  the  height  of  burst  curve 

with  zero  height  of  burst  is  well  defined.  For  absolute  heights  of  burst  less  than 
5/18 

1 m/kt  the  crater  volume  will  be  given  by 

V = V y5/610-0.49(H/Y5/18)  (23) 

o 


where  Vq  is  the  volume  of  a nuclear  crater  at  zero  height  of  burst  and  a yield 
of  1 kt,  and  Y is  in  kilotons.  By  contrast,  Figure  24  would  give 


V = V Y 10 
o 


-6.6(HOB/V1/3) 


which  is  valid  only  for  0<HOB/V1  < 0.02  which  includes  only  half  of  the  Pacific 
events  and  is  not  valid  for  any  negative  heights  of  burst. 

The  PACE  series  of  1000  lb.  TNT  events  conducted  at  the  Pacific  Proving 
Ground  permit  a direct  comparison  of  nuclear  and  high  explosive  height  of  burst 
effects.  Figure  26  shows  the  effect  of  height  of  burst  on  cratering  efficiency 

1 / O C /I Q 

for  both  V ' scaling  and  Y ' scaling.  The  curves  labeled  B and  C and  labeled 

E are  taken  from  Figures  24  and  25,  respectively,  with  extensions  at  depth  of 

burst  to  include  JANGLE  U.  Curves  A and  D,  respectively,  are  indicative  of 

the  slope  of  the  VA  and  Y°  height  of  burst  curves  for  high  explosive  as 

they  cross  zero  height  of  burst.  The  slopes  of  the  various  curves  at  zero  arc 

given  in  Table  8 along  with  their  intercepts.  We  note  that  the  slope  of  the 

5/18 

nuclear  curve  is  4.27  times  that  of  the  high  explosive  curve  for  Y scaling. 
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Figure  25.  Y5J8  scaling  applied  to  height  of  burst  effects  for  nuclear  events 


HOB 


Figure  2b.  Height  of  burst  curves  for  nuclear  and  high  explosive  events 


Table  8 

Relation  of  Height  of  Burst  to 

Efficiency 

Curve 

d (iog10  v/Yn) 

V/Yn  O 

HOB  = 0 

d (scaled  HOB) 

B 

- 6.60 

1.58 

x 104 

C 

-40.5 

5.62 

x 103 

E 

- 0.492 

1.15 

x 104 

A 

- 3.08 

1.09 

x 105 

D 

- 0.115 

2.79 

x 104 

0 


Hence, 


V = V Y5/6l<f°-115(H/Y  / \ (25) 

o 

We  are  now  in  a position  to  use  curve  E to  correct  for  height  of  burst 
effects  and  determine  the  effect  of  yield  on  crater  efficiency  for  Pacific  nuclear 
tests.  This  is  done  in  Figure  27.  It  can  be  seen  that  there  »s  considerable 
scatter,  but  the  least  squares  line  through  the  data  has  a slope  of  -.152,  very 
nearly  the  value  of- -.159  observed  in  Ottawa  sand.  We  can  also  determine  the 
source  effect  for  nuclear  sources  versus  TNT  sources  in  wet  coral  sand.  Extrap- 
olating to  1000  lb.,  we  see  that  TNT  is  2.75  times  as  effective  as  a nuclear 
device  at  zero  height  of  burst  (i.e.,  vq(S)  = 2.75).  This  is  40%  to  60%  lower 
than  the  value  obtained  by  comparing  craters  in  alluvium.  However,  in  view  of 
the  scatter  of  the  Pacific  data,  the  difference  may  not  be  significant. 

We  can  use  Schmidt's  tests  18-0  and  18-X  to  evaluate  height  of  burst 
effects  in  alluvium.  Based  on  these  two  tests  alone,  the  slope  similar  to  the 
values  listed  in  Table  8 of  a HOB-efficiency  curve  is  -0.133,  nearly  identical  to 
that  of  curve  D in  the  PACE  series.  This  agreement  leads  us  to  conclude  that 
the  two  tests  exhibit  a high  degree  of  internal  consistency.  The  reproducibility 
of  crater  volume  at  zero  height  of  burst  is  discussed  in  Section  4.4. 

4.4  PREDICTING  CRATER  VOLUMES 

The  relations  derived  in  the  last  two  sections  can  be  used  to  predict 

crater  volumes  to  be  expected  from  spherical  explosions.  As  an  example,  consider 

that  one  wishes  to  predict  the  crater  volume  from  a 1.25  gram  PETN  sphere  in 

alluvium  (density  1.57  Mg/m  ) with  a height  of  burst  of  -0.64  charge  radii  at 

345  g's.  (This  is,  of  course,  Schmidt's  18-0  test.)  Using  these  values,  plus  C-J 

data  from  Table  7,  we  use  Eq.  (21)  to  get  an  expected  volume  at  zero  height 

3 

of  burst  of  42.5 ± 15  cm  . Then  applying  a height  of  burst  correction  as  in  Eq. 

3 

(25)  we  get  an  expected  buried  volume  of  55  ± 20  cm  . 

This  does  not  compare  very  favorably  with  the  observed  volume  of  78 
3 

cm  . However,  the  scatter  between  shots  16  and  17  by  Schmidt  is  as  large  as 
the  scatter  among  many  alluvium  shots.  If,  however,  we  use  shot  16-0  or  16-X, 


3 

the  zero  height  of  burst  volume  would  be  about  54  cm  , giving  a buried  value 
3 

of  70  cm  , virtually  identical  with  that  observed.  This  illustrates  the  difficulty 
caused  by  variable  material  properties  in  cratering  experiments.  The  techniques 
used  to  prepare  the  test  bed  in  each  set  of  shots  - (16-0  + 16-X),  (17-0  + 17-X) 

3 

and  (18-0  + 18-X)  - were  slightly  different.  The  bulk  density  was  1.52  Mg/m 
fr  16-0  and  16-X,  1.61  Mg/m^  for  17-0  and  17-X  and  1.57  Mg/m^  for  18-0  and 
18-X. 


The  internal  consistency  of  the  two  events  18  can  be  checked  by  using 

Eq.  (25)  to  compare  their  observed  volumes.  The  difference  in  scaled  height  of 

5/18 

burst  of  the  two  events  was  1.28  m/kt  (18-0  was  at  0.64  charge  radii  and 
18-X  was  at  1.47  charge  radii,  both  below  the  surface.)  Hence,  we  would  expect 
that 


V18-X  = 1Q0. 115(1. 28)  = ]AQ 
V18-0 

3 3 

The  observed  volumes  were  78  cm  and  116  cm  in  a ratio  of  1.49,  only  6% 
higher  than  anticipated.  Hence,  we  conclude  that  the  tests  18  are  consistent 
with  other  height  of  burst  data  from  high  explosive  tests.  Crater  volumes  for 
each  set  of  tests  by  Schmidt  were  highly  reproducible  although  they  differed 
from  set  to  set  by  as  much  as  45%.  Collection  of  useful  data  in  centrifuge 
cratering  tests  will  depend  on  very  careful  attention  to  sample  preparation. 
Consistency  observed  within  a set  of  data  indicates  that  reproducible  materials 
can  be  obtained  for  a given  sample  preparation  technique. 

4.5.  OTHER  CRATER  MEASURES 

In  order  to  use  predicted  crater  values  for  targeting  or  vulnerability 

considerations,  it  is  necessary  to  predict  crater  radii  since  that  is  the  basis  of 

most  crater  related  vulnerability  measures.  Two  measures  of  crater  shape  are 

useful  in  this  regard  - the  aspect ' ratio,  R/d,  where  R is  the  apparent  crater 

2 

radius  and  d is  the  apparent  crater  depth,  and  the  shape  factor,  F = V/(ttr  d)  = 
3 

V(R/d)/iTR  . Our  discussion  will  involve  mainly  the  aspect  ratio,  although  a few 
thoughts  are  directed  toward  the  shape  factor. 
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A cursory  glance  at  data  for  R/d  from  Ottawa  sand,  wet  coral  sand, 

and  simulated  alluvium  with  high  explosive  and  alluvium  with  nuclear  explosives 

indicates  that  the  aspect  ratio  varies  markedly  from  one  material  to  another, 

and  also  with  height  of  burst.  If  the  aspect  ratio  is  plotted  against  scaled  height 

of  burst  for  a given  set  of  data,  a straight  line  can  be  fit  for  absolute  heights 

5/18 

of  burst  less  than  7 m/kt  . The  intercept  of  this  line  with  zero  height  of 

burst  is  then  taken  as  the  reference  aspect  ratio  for  that  data  set,  lR/d)x  for 

x ^ 

data  set  x.  Values  of  (R/d)Q  obtained  in  this  way  are  given  in  Table  9. 

3 x 

The  effect  of  yield  or  g Y on  (R/d)Q  is  negligible.  The  mean  of  7 high 

g tests  in  Ottawa  sand  at  zero  height  of  burst  was  4.50  which  is  only  11%  higher 

than  (R/d)^^*  listed  in  Table  9.  That  variation  is  less  than  one  standard 
o 

deviation  from  the  overall  reduced  aspect  ratio  versus  scaled  height  of  burst 
curve  derived  below. 

In  order  to  consider  R/d  data  from  several  environments,  we  introduce 
the  reduced  aspect  ratio  defined  as 


R 4. i\x  (R4joKE,a'- 

HOB  d'  HOB  (R/d)* 


Thus,  all  aspect  ratios  are  converted  to  values  consistent  with  nuclear  explosives 
in  alluvium.  When  reduced  aspect  ratio  is  plotted  against  scaled  height  of  burst, 
the  data  fall  into  a single  scatter  band.  Fitting  a straight  line  to  all  of  the 
data  at  once,  wc  find  that  there  is  a weak  dependence  of  (R/d)  on  height 
of  burst.  Specifically, 

(R/d)R0B  = 2.35  x 100-02'(HOB/Y5/18>  t 0.055  (27) 

where  the  indicated  uncertainly  is  one  standard  deviation.  This  line  is  shown  in 
Figure  28  along  with  the  data. 
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Table  9 

Zero  Height  of  Burst  Aspect  Ratio  of  Various  Crater  Data  Sets 


Date  Set 

Nuclear,  alluvium 
Boeing,  alluvium 
PACE,  coral  sand 
UDRI , Ottawa  sand 


Aspect  Ratio  (R/d) 
2.35 
2.60. 

2.77 

4.07 
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Figure  28.  Reduced  aspect  ratio  as  a function  of  scaled  height  of  burst. 


An  attempt  to  apply  Eq.  (27)  to  Pacific  nuclear  events  will  not  be 
successful  as  is  illustrated  in  Figure  29.  The  line  and  error  band  from  Figure  28 
are  included  along  with  data  from  JOHNIE  BOY  and  JANGLE  S.  Only  two  Pacific 
events,  CACTUS  and  LACROSS,  are  consistent  with  any  of  the  other  data.  Recall 
that  the  PACE  series,  which  were  1000  lb.  TNT  tests  at  the  Pacific  Proving 
Ground,  fit  Eq.  (27)  quite  well. 

The  variation  of  Pacific  nuclear  events  from  Eq.  (27)  is  primarily  a yield 
dependent  effect.  This  is  illustrated  in  Figure  30  which  is  a plot  of  the  ratio 
of  observed  (R/d)  to  that  calculated  by  Eq.  (27).  A straight  line  has  been 
fit  to  the  data  with  a standard  deviation  as  indicated  by  the  dashed  lines.  (The 
tank  shots  were  ignored  in  the  fit.) 

The  deviation  of  aspect  ratios  of  Pacific  events  at  higher  yields  is  not 

(12) 

a new  result.  Numerous  explanations  such  as  water  washing  , crater  slump- 
ing^41 \ and  dynamic  liquification^  ^ have  been  offered.  Any  explanation  should 
be  consistent  with  the  observed  dependence, 

<R/d>obs/(R/d)pred  “ Y''15  I0-5-  (28) 

A completely  empirical  observation  is  noted  which  may,  eventually,  lead 
to  a resolution  of  the  wide  scatter  observed  in  Pacific  nuclear  crater  dimensions. 
Events  whose  reduced  aspect  ratio  is  higher  (lower)  than  the  line  in  Figure  30 
have  larger  (smaller)  volumes  than  predicted  from  the  line  in  Figure  27.  This 
is  clearly  indicated  in  Figure  31.  Since  the  deviations  from  Eq.  (28)  are  strongly 
correlated  with  those  from  volume  data,  we  must  conclude  that  anomalies  in 
shape  and  volume  are  due  to  a single  cause. 

If  the  shape  factor  F is  plotted  against  scaled  yield  for  the  Ottawa  sand 
tests  shown  in  Table  6,  a trend  to  higher  F at  higher  scaled  yield  can  be 
ascertained.  Numerically,  it  is  found  that  F = (0.58+  0.03)  g^W/(U^i ) 

Because  of  this  slight  yield  dependence  of  F,  the  crater  dimensions  R and  d do 
not  quite  scale  as  V1^.  Rather,  we  have  (R,d)  a _ ^1/3.6. 

One  final  point  on  crater  shape  factors  (F)  will  be  made.  There  is  little 


72 


Figure  29.  Reduced  aspect  ratio-scaled  heioht  of  burst  effects  for  nuclear  tests. 
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Y(kt) 

Figure  30.  Yield  dependence  of  deviation  of  R/d  for  Pacific  events. 


volumes  and  aspect  ratios  from  mean 
values  taken  from  Fioures  27  and  30. 


75 


height  of  burst  dependence  of  F seen  in  the  PACE  series,  and  values  are  similar 
to  those  of  lower  yield  Pacific  nuclear  events,  as  well  as  JOHNIE  BOY  and 
JANGLE  S (all  are  0.5±  0.1  except  SEMINOLE  (a  tank  shot))  which  has  F = 0.67. 
However,  events  with  yield  over  100  kt  have  F = 0.26  ± 0.04,  except  for  ZUNI 
and  TEWA.  ZUNI  and  TEWA  exhibit  the  largest  volume  and  aspect  ratio  anomaly 
(of  the  non-tank  events)  as  well  as  being  anomalous  in  F.  However,  a correlation 
similar  to  that  shown  in  Figure  31  is  not  apparent. 
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tests  shown  in  Table  6,  a trend  to  higher  h'  at  mgner  scaieu  yitiu  can  i/t 
ascertained.  Numerically,  it  is  found  that  F = (0.58 1 0.03)  g3W/(U66 ) 0,0045 
Because  of  this  slight  yield  dependence  of  F,  the  crater  dimensions  R and  d do 
not  quite  scale  as  V^3.  Rather,  we  have  (R,d)  a W3,33®  = W^3’®. 

One  final  point  on  crater  shape  factors  (F)  will  be  made.  There  is  little 


Figure  31.  Deviation  of  observed  Pacific  crater 
volumes  and  aspect  ratios  from  mean 
values  taken  from  Figures  27  and  30. 


